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ABSTRACT. The total energy in 3D MHD turbulence is shown to cascade over an inertial range the 
lower limit of which is a modified Taylor length scale. Direct cascades are also investigated for 
energies transported by distinct mechanisms which together comprise the total energy. These include 
the convection of total fluid energy by the velocity field, the kinetic to kinetic energy transfer advected 
by the velocity field, and the magnetic to magnetic energy transfer advected by the velocity field. Flux 
locality for the cascading quantities is also addressed. A scenario in which the inter-medium energy 
transfer is predominantly from kinetic to magnetic and a note regarding the cascade of kinetic energy 
for non-decaying fluid turbulence are included. 

1. INTRODUCTION. 

Turbulence remains an enigmatic feature of modern fluid and plasma dynamics. In the case of 
3D incompressible fluids, the first positive rigorous result pertaining to the existence of energy 
cascades was achieved in Fourier space by Foias, Manley Rosa, and Temam IfTTll . A pioneering 
step was recently taken in [ToU where a mathematical apparatus was developed to study turbu- 
lence working entirely in physical scales. In this setting it was shown that ensemble averages of 
localized space-time averages (localized to certain coverings of an integral domain) of the energy 
flux through shells are positive and nearly constant. Consequently, the transport of energy is, in 
an averaged sense, directed from larger to smaller scale structures. Interestingly, locality of the flux 
in this context is naturally affirmed. More precisely, the locality is derived dynamically as a direct 
consequence of existence of the turbulent cascade in view, featuring comparable upper and lower 
bounds throughout the inertial range; in contrast, the previous locality results were essentially 
localized kinematic upper bounds on the flux, the corresponding lower bounds being consistent 
with turbulent properties of the flow |fT6l l8ll. The formal similarities between the magnetohydro- 
dynamical system and the Navier-Stokes equations allow the theory of the former to be informed 
by that of the latter. The present work proceeds in this spirit. Here, the methodology of Dascal- 
iuc and Grujic is applied to establish the existence and locality of various energy-level cascades 
in turbulent MHD flows over suitable inertial ranges. An additional comment is included re- 
garding preferential inter-field energy transport. A complementary work by the authors |6| is in 
preparation and will deal with kinetic and magnetic enstrophy cascades in 3D incompressible MHD. 
The two works combined will constitute a solid foundation to the theory of 3D MHD turbulent 
cascades/ transport in physical scales. 

Magnetohydrodynamical turbulence research has been energized by recent solar wind observa- 
tions as well as the ability to generate more complex and more accurate numerical models. The 
picture painted by these advances is, however, far from conclusive. The most recent theories indi- 
cate energy cascades in certain settings and posit power-laws with dependencies on the viscosity 
and resistivity, as well as global quantities associated with the flow fl3j |22j 01 |2lJ [19| . The first 
purpose of this paper is to rigorously establish the cascade of total energy by studying the total 
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energy flux in a suitably statistical manner across a range of scales - the inertial range - for which 
the lower bound is a corrected Taylor micro-scale, the correction of which is dependent on param- 
eters arising from the employed mathematical apparatus as well as the magnetic Prandtl number. 
This result follows directly from the velocity-magnetic field formulation of the 3D incompressible 
magnetohydrodynamic equations and is physically reasonable for all forms of MHD turbulence 
where the magnetic Prandtl number is not significantly smaller than one (i.e. rj < v). In particular, 
it applies in astronomical settings such as the Solar wind and the interstellar medium where n is 
often so small that the non-resistive limit is commonly used as an approximation 0. Interest- 
ingly in contrast to current numerical and phenomonological theories, no appeal is made to the 
existence of a strong magnetic guide field and we thus assert that even in the absence of a strong 
background magnetic field, 3D MHD turbulence exhibits an inertial range over which a cascade 
of total energy transpires. 

Due to the coupling between the magnetic and kinetic fields, each of which is imbued with its own 
energy there are a number of transfer mechanisms by which energy can 'flow' between scales. In 
particular, energy can remain tied to the kinetic field or the magnetic field but it can also be trans- 
ferred between the two fields-i.e. from the kinetic field to the magnetic field or vice versa. The 
secondary purpose of this work is to identify conditions under which the distinct energy transfers 
are, in an ensemble averaged sense, directed from larger to smaller scales. Particular attention is 
paid to the term responsible for the transfer of energy locked within the velocity field. This is the 
total fluid energy flux. Interestingly, the numerically and observationally supported phenomenon 
of dynamic alignment - the scale dependent depletion of the angle between the direction of the ve- 
locity and magnetic fields |3H22j|4j|23 _ can De interpreted to obtain a (weak) geometric depletion 
of the non-linearity complementary to the fluid energy flux. We investigate the role this plays in 
the relevant cascade. 

In |Tl3ll , the ensemble averaging mechanism which enables our study of turbulence in physical 
scales was interpreted to study non-flux type quantities (the vortex stretching term) and establish a 
mathematical evidence for the persistence of the lengths of vortex-filaments across a (longitudinal) 
range of scales. A similar argument can be applied to establish a scenario in which the inter- 
field energy exchange is predominantly oriented from the velocity field to the magnetic field. To 
achieve this we provide a dynamic estimate which entails the stretching effect of the velocity field 
on the magnetic field lines is predominantly one of elongation (the magnetic field is energized at 
the expense of the velocity field) as opposed to diminution (in which case the velocity field would 
be energized). 

It is important to mention that pursuing a rigorous study of turbulent cascades in physical scales is 
a break from the traditional approach (9j[l8|. Despite the dependence on properties of turbulence 
which are emergent in physical space present in the Kolmogorov dimensional analysis, the predom- 
inant setting for the study of these cascades has been in Fourier space, i.e., in the wavenumbers; 
this required postulating (in a statistical sense) isotropy and homogeneity of the flow lTl7l|T6l l8ll. 
A direct physical interpretation of scale as a platform to understand cascades has traditionally been 
rejected. The crucial problem for a physical interpretation is that the cascade does not occur lo- 
cally - the energy flux cannot be expected to be, in an averaged sense, inwardly oriented on any 
particular shell. In the present framework, this is overcome by an ensemble averaging process 
whereby any local detractions from the essential inward-ness of the energy flux are shown to be 
insignificant and thus, in a statistically significant sense, the net 'flow' of energy is from larger to 
smaller scales. 

Several subjects treated in the present manuscript are noteworthy in the context of the existing 
literature regarding fluid and plasma turbulence in physical scales. First is the development of 
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cascades for fluxes which exclude the pressure flux. Because analysis of the pressure flux is fun- 
damentally different than other fluxes appearing in both MHD and NSE, additional steps must be 
carried out to interpolate ensemble averages of localized quantities between integral scale quan- 
tities. Our present treatment of the fluid pressure flux provides an avenue for the study of non- 
decaying turbulence which will be investigated in a future work. Second, the relationship between 
cascades for the dimensional and dimensionless forms of the equations is investigated and found 
to be equivalent. This is technically beneficial because, by carrying out analysis after passing to 
the dimensionless form, issues of dimensionality are tidily avoided. 



The main purpose of this section is to describe how ensemble averaging with respect to (K~i,K 2 )- 
covers of an integral domain B(0, Rq) can be used to establish essential positivity of an a priori sign- 
varying density over a range of physical scales associated with the integral domain (cf. IflOl ITT1 
IT4l [El). The application to turbulence is establishing the positivity of certain inward directed 
flux densities - i.e. the cascade is uni-directional from larger to smaller scales - as well as the 
near-constancy of the averaged densities - i.e. the space-time averages over cover elements are all 
mutually comparable - across a range of scales. 

The ensemble averages will be taken over collections of spatio-temporal averages of physical den- 
sities localized to cover elements of a particular type of covering - a so called (K\, K 2 )-cover - 
where the cover is over the region of turbulent activity. For simplicity, this region will be taken as 
a ball of radius Rq centered at the origin and, to reflect the turbulence literature, is henceforth re- 
ferred to as the integral domain (also known as the macro-scale domain). The time interval on which 
we localize is motivated by the physical theories of turbulence and assumed to satisfy, 



The {K\, K 2 ) -covers are now defined. 

Definition 1. Let K\, K 2 € N and < R < Rq. The cover of the integral domain B(0, Rq) by the n balls, 
{B(xi, R)}i=i is a (Ki, K 2 )-cover at scale R if 



and, for any x e B(0, Rq), x is contained in at most K 2 balls from the cover. 

In the hereafter, all covers are understood to be (Ki, K 2 )- covers at scale R. The positive integers 
K\ and K 2 represent the maximum allowed global and local multiplicities, respectively. 

In order to localize a physical density to a cover element we incorporate certain refined cut-off 
functions. For the cover element around the point X{, let 4>i(x, t) = n(t)tp(x) where r\ € C°°(0, T) 
and ip € C^(B(xi,2R)) satisfy 



2. (Ki,K 2 )-Covers and Ensemble Averages. 





< n < 1 



t] = Oon (0,T/3), 



t? = 1 on (2T/3, T) 



\dtv\ ^ C 



(1) 



and 



o < i> < i 



ip = 1 on B(xi, R) 



I^VI ^ Co 



ipp ~ R 



\difW\ Cp 

^2p-l - R 2 ■ 



(2) 



where 3/4 < 5, p < 1. 
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By 4>q we denote the cut-off function associated with the integral domain - the ball centered at 
x = of radius Ro - satisfying the above properties. 

Comparisons will be necessary between averaged quantities localized to cover elements at some 
scale R and averaged quantities at the integral scale. To accommodate this we impose several 
additional conditions for Xj near the boundary of B(0,R ). If B(xi,R) C B(0,Rq) we assume 
ip < ip - Alternatively, when B(xi, R) <£ B(0, R ) we stipulate that ip = 1 on B(x , R) n B(0, R ), 
satisfies (O, and we additionally have: 

ij) = ipo on the intersection of S(xq, Ro, 2Rq) and the cone with apex at the origin 
and with boundaries passing through the intersection of the circle centered at the 
origin of radius Rq and the boundary of B(xi, R), 



The above apparatus is employed to study properties of a physical density at a physical scale R 
associated with the integral domain B(0,Rq) in a manner which we now illustrate. Let 9 be a 
physical density (e.g. a flux density) and define its localized spatio-temporal average on a cover 
element at scale R around 



1 JO n JB(xi,2R) 

where < 5 < 1, and let (@}r denote the ensemble average over localized averages associated 
with cover elements, 



Examining the values obtained by ensemble averaging the averages associated to a variety of 
covers at a fixed scale allows us to draw conclusions about the flux density 9 at comparable and 
greater scales. For instance, stability (i.e. near constancy) of {(6).r} indicates that the sign of 9 
is essentially uniform at scales comparable to or greater than R. On the other hand, if the sign 
were not essentially uniform at scale R, particular covers could be arranged to enhance negative 
and positive regions and thus give a wide range of sign varying values in {(6)^}. In order, then, 
to show the essential positivity of an a priori sign varying density 9 at a scale R, it is sufficient to 
show the positivity of all elements of {(0},r}. 

Conversely, if 9 is an a priori non-negative density, then ensemble averages are all comparable 
to the integral scale average across the range < R < Ro. We make this notion precise in the 
following lemma. 

Lemma 2. Let f(x,t) € Lj oc ((0,T) x ]R 3 ) be non-negative. Let {xi}f =1 be centers of elements of a 
(Ki,K2)-cover of B(xo,Rq) at scale R < Rq. Setting 



and, 



i/j = on the intersection of the three sets B(0, R ) \ B(x h 2R), 5(0, Rq, 2R ), and 
the outside of the cone with apex at the origin and boundaries passing through the 
intersection of the circle centered at the origin of radius Rq and the boundary of 

B{x u 2R). 






i=l 








and 
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, F < (F) R < K 2 F . (3) 



Proof. Recalling that 4> Xi ,R < 4>o an d the definition of (K±, ^"2) -covers we have that, 

(F) R = ^^ -ij y f{x,t)Yj Xi)R {x,t)dxdt 

~tJ 1$JP J f( x ^) K Mx,t) dx dt = K 2 F , 



and, 



(f) r = i y _l y /(x, t) £ ^^(x, *) dt 

1 /" T 1 1 R 3 f 1 
-fj WlT^l fMM*,t)dxdt = -F . 



□ 



For additional discussion of {K\, K 2 ) -covers and ensemble averages, including some computa- 
tional illustrations of the process, see (14] . 



3. 3D Incompressible MHD Equations. 

Our mathematical setting is that of weak solutions to the 3D magnetohydrodynamic equations over 
M 3 (cf. EH for the foundational theory). Define V = {/ G L 2 (M 3 ) : V-/ = 0} (where the divergence 
free condition is in the sense of distributions) and let V be the closure of V under the norm of the 
Sobolev space, (// 1 (M 3 )) 3 , and, H, the closure of V under the I? norm. By a solution to 3D MHD 
we mean a weak (distributional) solution to the following coupled system (3D MHD): 

u t - M11 + (u ■ V)u - (b ■ V)6 + V(p + \b\ 2 /2) = 0, 
bt - r]Ab + {u ■ V)6 - {b • V)tt = 0, 
V • u = V • b = 0, 
0) = 1*0(3;) G V, 
b(x,0) = b (x) G V, 

where r\ and ^ are the magnetic diffusivity and kinematic viscosity respectively and p(x, t) is the 
fluid pressure. 



Our present work utilizes suitable weak solutions for MHD, our definition of which follows that in 

ED]. 
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Definition 3. The functions u, b, p constitute a suitable weak solution to 3D MHD if they are a weak 
solution and further satisfy the following generalized energy inequality: 

J J (v\\7u(x,t)\ 2 + rj\\7b(x,t)\ 2 )(f>(x,t) dx dt (4) 
< \J J (\u{x,t)\ 2 + \b{x,t)\ 2 )fa{x,t)dxdt 

+ \ I j ( u \ u ( x,t }\ 2 + r l\ b ( x ^ t )\ 2 ) A< f ) ( x ' t ) dxdt 

+ \j j ^ x,t ^ 2 + H X M 2 + 2 v{x,t)){u{x,t) - V(j){x,t)) dxdt 

- J J (u(x, t) ■ b(x, t))(b(x, t) ■ V(p(x, t)) dx dt, 
for a.e. T G (0, oo) and any non-negative 4> G Co°(M 3 x [0, oo)). 

Existence of suitable weak solutions for MHD is proven in [20] using an adaptation of the tradi- 
tional method for NSE first presented in 0. We will make use of certain regularity properties 
which are evident from the existence theorem given in |20l . 

Proposition 4. For uq, bo € H and uq e W 4 ^ 5 ' 5 ^ 3 , suppose (u, b,p) constitutes a suitable weak solution 
to 3D MHD. Then, (u, b,p) satisfies, 

u,b£ L°°(0,oo;F), u,b G L 3 (0, oo; L 3 (R 3 )), 

Vu,V6 G L 2 (0, oo; L 2 (M 3 )), p G L 3/2 (0, oo; L 3/2 (M 3 )). 

The fact that suitable weak solutions only satisfy a generalized energy inequality (as opposed to 
equality) introduces the possibility that energy is dissipated not only by viscosity or resistivity 
but also by singularities. In the case that the weak solution in question is regidar, equality is 
attained in the generalized energy inequality ((U, and the possibility of potential loss of flux due 
to singularities is eliminated. To streamline discussion we establish cascades for regular solutions 
and include a result for the non-regular case only in the context of the cascade for the modified 
(due to energy loss from possible singularities) total energy flux (c.f. Section @]). The study of 
the possible energy loss due to singularities is itself an interesting subject but the case of MHD is 
not sufficiently distinct from that of NSE (which can be found in |H0"|0 to justify an independent 
exposition. 



4. Total Energy Cascade. 



The total energy flux through the boundary of the ball B, over the interval (0, T), is given (cf. (3j) 

-/ / (\u\ 2 + \b\ 2 + 2p) h ■ u dx dt - I [ (u ■ b)(h ■ b) dx dt, 

2 Jo Job Jo Job 

where n is the unit normal vector directed inward. Our analytic results are enabled by substitut- 
ing the inwardly directed vector field V</> for h where is a refined cut-off function for the ball 
B(xq, R). Our localized, space-time averaged total energy flux into the ball centered at the spatial 
point xq of radius R, over the interval (0, T), is then defined as, 

F x ,R-=\J Q J(\ u \ 2 + \ b \ 2 + 2p)(u-V(f>)dxdt- J J(u-b)(b-V(f>)dxdt. 
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The following remark on the genesis of the last term, the advection of the cross-helicity via the 
magnetic field, is informative. The transfer of magnetic to kinetic energy is driven by the Lorentz 
force while the stretching of the magnetic field lines is responsible for the transfer of kinetic energy 
to magnetic energy. Since these are complimentary, the sum cancels in the non-localized case. 
Locally, however, we obtain a flux-type term, 

((b ■ V)b • <j)u + (b ■ V)u ■ <j)b) dx dt = - J(u-b)(b-V(p)dxdt. (5) 

The appearance of this term is interesting because it is only non-zero if there is some degree of 
non-locality in the energy transfer between the two fields. 

In the following we work in a fixed integral domain, B(xq, Ro), with associated cut-off (f>o = 4>x Q ,R Q - 
Certain integral domain quantities will be used to determine lower bounds on the inertial ranges 
over which our cascades are rigorously shown to persist. Our labeling scheme follows. The inte- 
gral scale space-time averaged kinetic and magnetic energies are defined in terms of a technical 
parameter, 5, as, 

e u = e u (8) = ± £ j^j l -\u\mdxdt and e b = e b (5) = 1 £ ^ J ' \\b\ 2 <j>l dx dt, 

and the integral scale time-space averaged enstrophies by, 

E% = v^J^ J \Vu\ 2 (f>o dx dt and E b = r,^ J |V6|^ dx dt. 

Regarding notation, we indicate the combined kinetic and magnetic energies or enstrophies by 
omitting the superscript (i.e. eo := e^ + e b Q ) and, because eo is decreasing with 8, take liberties 
suppressing the dependence of eo on 5 with the understanding that the indicated quantity is that 
associated with the smallest appropriate value for 5. 

We will establish that the cascade persists over a range of scales bound above by the integral scale 
and below by a modified Taylor micro-scale. The Taylor micro-scale is (c.f. l3ll), 

/ \ 1/2 

The modification will depend in part on the magnetic Prandtl number, denoted by Pr, a non- 
dimensional number given by the ratio of kinematic viscosity to magnetic diffusivity, i.e., 

_ v 
Pr = -. 

V 

Essentially its roll in the modification is to incorporate information about the magnetic diffusivity 
which is absent from our prescribed time scale (recall T > R\j v) and the Taylor micro-scale. 

Theorem 5. Assume u and b are suitable weak solutions to 3D MHD with sufficient regularity that equal- 
ity holds in Let {xi}f =1 C B(xq, Rq) be the centers of a {K\,K-})-cover at scale R where the cut-off 
functions are defined with T > R 2 jv. For a scale and cover independent positive parameter 

1/2 

0:-- 



2CK 1 K 2 {1 + Pr- 1 ) 

where C is a constant determined by structural properties of 3D MHD and our cut-off functions, if 

t < f3Ro, 

then, 

1 -E <(F E ) R <2K 2 E , 



2K - r ir 
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provided R is contained in the interval [t//3, Rq]. 

Proof. Assuming the premises above and in virtue of 0$, we have for any cover element that, 

F% iR > £ J {v\Vu\ 2 + 77| Vb\ 2 )^, R dx dt 

-\J j\\u\ 2 + \b\ 2 )d t ^ XiyR dxdt + ^ J J{v\u\ 2 + r ] \b\ 2 )A ( f )Xi , R dxdt . 
Recalling the property of our cut-off function, 



as well as the fact, 



i a j. i ^ ^xi,R 
\Ot<Pxi,R\ < CO^-, 



-< — = Pr^- 
T ~ R 2 R 2 ' 



we conclude that, 

^ lo l (|U ' 2 + ^ dt<t>Xi ' R ^ ^ ~ C °^2 £ j \ U \ 2 € % ,R dx dt + C Pr ^2 f / \ b K it R dx dL 

Our cut-off functions also satisfied, 
and, consequently, 

\ J J ("H 2 + vlb^A^x^R dx dt < J J \ u \ 2< Px P ~R dx dt + °o^2 J J l^l 2 ^,/? 1 dx dt - 
Noting that p > 2p — 1 we have, upon combining the above estimates, that, 

F% iR > £ j (u\Vu\ 2 + n\ Vb\ 2 )^ R dxdt-^£ J Hu| 2 + r?(l + Pr)\b\ 2 )^J dx dt. 
Using Lemma we observe that, 

K 2 E >(^J q ^3 J (^|Vu| 2 + ??|V6| 2 )^, !ii? dx dtj > ^-Eq, 



and, 

T 



ff Q Jp J\^H 2 + V(l + Pr)\b\ 2 )^dxdt^ R <uK 2 e u 0+V (l + Pr)K 2 e b 

<vK 2 {l + Pr~ l )e Q . 



We can thus interpolate the ensemble average between integral scale quantities as, 

1 „ c K 2 {l + Pr~ l ) E cpK^l + Pr- 1 ) 
—E - v — e <(F ) R < K 2 E + v — 2 e . 

Note that the upper bound follows in virtue of the assumed regularity (i.e. equality in ((U)) and is 
the only place where this assumption is used. In particular, the lower bound holds for non-regular 
solutions. 

Continuing, we now specify a value for (3, the modification to the inertial range, to be, 

/ i \ 1/2 



2c Q K 1 K 2 (l + Pr- 
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Using the inertial range condition that t/ (3 < R < Rq, we have, 

cog 2 (l + Pr- 1 ) 1 

" # 60 * 2i^°" 

The final lower bound for the ensemble average then becomes, 



',F E )„>^E . 



R ~ 2K X 

The upper bound follows trivially with our definition of ft and we conclude that, 

— !— E < (F E ) D < 2K 2 E . 

□ 

Remark 6. T/ze dependence of the length of the inertial range on Pr has consequences for when the above 
result is most physically relevant. The correction to the Taylor micro-scale is minimized when Pr^ 1 < 1, 
that is, when r/ < v. This is exactly the scenario for turbulence of plasmas in astronomical settings such as 
the solar wind and the interstellar medium. 

We briefly discuss how to obtain a result for possibly non-regular suitable weak solutions. As 
noted in the proof, the assumption of regularity was only used in establishing the upper bound on 
the ensemble average. This is the issue we now address. The physical cause of a strict inequality in 
the generalized energy inequality is interpreted as the loss of energy due to possible singularities. 
This will be denoted by F?° (or F£° R if is the cut-off for a {K\ , i^-cover element) and is defined 
as the value that fills in the inequality (H)), 

J j\v\Vu\ 2 + T]\Vb\ 2 )<p dx dt = ^ J j \\u\ 2 + \b\ 2 )<p t dx dt + i J J \v\u\ 2 + r]\b\ 2 ) dx dt 

° ° (6) 

+ \ J J(\ u \ 2 + 2 \ b \ 2 + 2 P)( u ' V( l ) ) dxdt - J q j {u-b){b-V(j>) dxdt- F^. 

To account for the possible strictness of the above inequality we replace the fluxes considered 
previously with fluxes modified to include the possible loss of energy due to singularities. For 
example, in the case of the total energy flux, we establish an interpolative bound on ensemble 
averages corresponding to the localized modified total energy fluxes, 

rpE,oo _ r?E _ poo 
r 4> ~ r 4> r 4>- 

Positivity and near-constancy results are then given in terms of the modified flux specified above. 
The cascade of total energy modified due to possible non-regularity is then given by the following theo- 
rem, the proof of which is identical modulo a substitution of F^ '°° for F E to the proof of ©. 

Theorem 7. Assume u and b are suitable weak solutions to 3D MHD. Let {xj}" =1 C B(xq, Rq) be the 
centers of a (Ki, K-i)-cover at scale R where the cut-off functions are defined with T > R 2 jv. For a scale- 
and cover-independent, positive parameter, 

1/2 

0:-- 



2CK 1 K 2 {l + Pr- l ) i 

where C is a constant determined by structural properties of 3D MHD and our cut-off functions, if 

t < f3R , 

then, 

1 -E <(F E >°°) n < 2K 2 E , 



2Ki u " x /R 
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provided R is contained in the interval [t//3, Rq\. 



5. Total Fluid Energy Cascade. 

Our attention is now turned to establishing conditions under which a cascade of the total fluid 
energy will exist. The localized flux quantity of interest is, 



/(H 2 + 2p)(u-V0)dxcft = -^ J((u-V)u + p) 



u) dx dt. 



Results are ultimately intended for the dimensional form of the equations (as was the case for 
the cascade of total energy). Our estimates, however, will substantially benefit from switching 
to the dimensionless form of the equations. Some general remarks illustrating this process are now 
included. Starting with a dimensional problem where the integral domain has radius Rq (the char- 
acteristic length scale) and v is the kinematic viscosity (these will be our fundamental dimensions), 
we will establish the existence of a dimensionless cascade over an inertial range given by, 

R ) Eq P ' 

where the dimensionless analogues to the dimensional solutions and variables are indicated with 
asterisks and the analogous normalized integral domain quantities are, 



1 



(\u*\ 2 + \b*\ z )<j> dx*dt 



*\2\ 



E* 



and f3 is a dimensionless constant. 



-(|Vu*| 2 



W| 2 )0 O dx*dt*, 



This has consequences for the dimensional setting. In particular, we have the quantitative rela- 
tionship, 



— — Ur,——, 



Eo 



^0 



which follows from a change of variable and the chain rule. Furthermore, setting our characteristic 
time scale to be T = R^u , we see the integral scale kinetic energy relates to its dimensionless 
analogue in the following way, 



e 



TR 3 n 



T 



o Jo 



\u\ cf>Q dx dt 



TRl 



o Jo 



R 2 , u(x,t) 

T 2 ~u7~ 



4>o(x, t) dx dt 



TRl 



o Jo 

T 



^\u*(x*,t*)\ 2 Mx*,t*) R* dx* T dt* 



R l 

f2 



3. P u* 

T 2 e 







R\ 

u*{x\t*)\ 2 (j) {x*,t*) dx* dt* 



Similarly, 



T 2 e J 



E - jQ. E i and ' E - J^ E 
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For the flux presently of interest, total (fluid) energy, setting, 

F Xt , R = ^3 J J((u ■ V)u + Vp) • (<pXi,Ru) dx dt, 

K,r = || J J ((«*(**, O ■ V.)u*(x*,t*) + V*p*(x*,t*)) ■ <t> XuR (x*,t*)u*(x*,t*) dx* dt*, 

we have, upon transforming the former into the latter, that 

p _ p* 

r Xi,R — rp 2 r Xi,R- 

Assuming now the existence of the cascade for the dimensionless system we can recover a cascade 
for the dimensional scenario as we now illustrate. We will see theorems pertaining to dimension- 
less cascades which assert that, for certain dimensionless quantities f3 u and fa, if, 

p / e «* \ V 4 p / J>* \ 1/4 



<R<R and < R < R , 

then, 

^E* < (F*)r < 2K*E* . 

By the quantitative relations identified above, the consequence for the dimensional scenario is, 
if, 

<R<R and -g- M- < R < Rq, 



then, 

^E < (F) R < 2vK*E . 

One last comment regarding the interplay between dimensional and dimensionless structures is 
that, upon making a change of variables with x* = x/Rq and t* = t/T and applying the chain 
rule, our refined cut-off functions are seen to enjoy the following bounds, 

|V^(x*,t*)|<^co|^*,t*)|" and \<f> t * (x* , t*)\ < c |^(x*, t*)\. (7) 

The above discussion justifies a switch to the dimensionless setting which we presently make. We 
subsequently suppress the asterisks used above to indicate non-dimensionality noting that for the 
remainder of this section all functions and variables are dimensionless. 



Following ||23|, dimensionless 3D MHD is, 

d t u--?-Au = -(u-V)u-Vp-s( V— -(b-V)b 
Re \ 2 

d t b- — Ab = V x (b x u), 
Rm 

V • u = V • b = 0, 

u(x,0) = u (x) G L 2 (M 3 ), 
b{x,0) = b (x) G L 2 (R 3 ), 

where Re and Rm are the Reynolds and magnetic Reynolds numbers respectively and the non- 
dimensional number, S, is defined in terms of the Hartmann number, M, to be S = M 2 / (ReRm). 

In order to obtain a local cancellation between coupled non-linear terms, we assume sufficient 
regularity so that the following derivation is justified. Taking the scalar product of the equation 
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of motion by <fiu and the induction equation by S4>b, integrating, and rearranging we obtain a 
flux density for total (fluid) energy into the ball B(xi/Ro, R/Rq) with associated cut-off function 

J / Q' U ' 2+P ) ' ^ V ^ ' u ) dx dt = Jfef Q J l Vu | 2( ^ dx dt + J J \^b\ 2 <pdxdt 



- S J J V(f> ■ (b x (u x b)) dx dt 

~ t: / / \ u \ 2 4>t dx dt — — — [ [ \u\ 2 A(f> dx dt 
2 Jo J 2i?e J J 

Bounds for the lower order terms on the right hand side follow. The last four terms are bounded 
in a manner similar to that seen in the proof of 10. Here, however, we acknowledge the change of 
variable and cite the properties of our refined cut-off functions (c.f. 10) as well as the chain rule to 
obtain, 

^ Jo I ' U ' 2 ^ dX ^ ~ C ° Re Jte (ir) lo I l n ' 2 ^ 4p ~ 3 dx dt > (8) 

1 P 1 /D.\2 yl 



±-fJ\ U fA,dxdt<col^) I f\u\ 2 ^dxdt, 
| £ J \bU t dx dt < CoRm^- f^j jf J |&|V dx dt, 
jf 1 J |6| 2 A</> dx dt < c -^ (j^j 2 jf * f \b\ 2 ^ dx dt. 



2Rm 

Repeatedly using Holder's inequality, the Gagliardo-Nirenberg inequality, and Young's inequal- 
ity, the following bound on the remaining term is evident, 



CqS ^ f 1 [ cjf{b x(ux b)) dx dt < c S^ f 1 I (|n| 1 / 2 ^- 3 / 4 )(|6| 1 / 2 )(|u| 1 /2| 6 |3/2 (/) 3/4 ) dx dt 



to 

R Jo J ~ " R j 



< c S^ f 1 \\u<?o-^\\y 2 \\b\\\' 2 ||(|u|V>|3/y/> dt 



R 







Rpf M^Rm 1 / 2 ^/ M 3 / 2 
- C ° R V Re 1 /* J V^e 3 / 4 ^ 3 / 2 

1 |n^- 3 / 2 ||f||6||f||n||f||V(6^/ 2) || 3 / 2 , t 

M 2 Rm 2 R% , rl 
t t Jo 







<&> J=£ a (3U P ||«||l)(sup||t||l) / |K 4 '- 3 ||l<ii 







< Co (A/iJm) 2 ( sup sup ||b|||) -i--|a r dt 



+ ARm 
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Taking ensemble averages and passing where appropriate (via Lemma 12)) to the quantities eft, e h Q 
En, and Ek, we have , 




where in the last line we have set, 



C u = cq\1 + Re + (Afiim) 2 (sup|H|2Sup ||6|||)J and C b = c (2 + Rm) 
Similarly, an upper bound is, 



(F)n < K 2 ±- e E% + C„K 2 -L§eZ + + CUT, Ag^. 



Defining now the parameters for a correction to the inertial range by, 




1/4 



1 



) 



1/4 



and /3f, 



2K 1 K 2 C b 



we have justified the following theorem. 

Theorem 8. Let {xj}™ =1 c B(xq, Rq) be the centers of a (Ki,K 2 )-cover at scale R. For j3 u and j3 h defined 
above, if 



Remark 9. It is worth noting that the condition triggering the cascade is essentially a requirement that the 
gradients of the velocity and the magnetic fields are large (averaged, over the integral domain) with respect 
to the fields themselves; this will hold in the regions of high spatial complexity of the flow (the correction 
parameter j3 depends on certain a priori bounded quantities; however, none of them involve gradients). 

Remark 10. A particular scenario where the above result is pertinent is one where S/Rm ~ 1 and the 
integral scale magnetic energy dominates the integral scale kinetic energy. The effect of the above assumption 
is to free the parameter j3 u from its dependence on Re by substituting for the bound the bound 



This allows for very large values of Re and is applicable to settings involving the confinement of a liquid 
metal by a strong magnetic guide field |24| . Without these assumptions on the flow the physical relevance 
of the above result is diminished by the fact that the dependences of j3 u and f3 b on the fluid and magnetic 
Reynold's numbers result in an inertial range which decreases in length as Rm and Re increase. 




then for scales R where t//3 < R < R ,we have, 
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An Improvement via Dynamic Alignment. Interestingly, the possible detraction from the posi- 
tivity of the ensemble average by the remaining flux terms from the generalized energy inequal- 

ity, 

N Xi>R := J J \b\ 2 (u ■ V<t> Xi , R ) dxdt- J j\u- b)(b ■ V<t> Xi , R ) dx dt, 

is (partially) balanced by an observationally and numerically supported geometric depletion in the 
form of scale dependent dynamic alignment (c.f. in particular Q as well as |5H22H2TJ|2]). Essentially, a 
decreasing length scale corresponds to an increase in the alignment or anti-alignment of the angle 
between u and b. It is a matter of vector calculus to rewrite N Xil R as, 

N XitR = J J V(f> Xt>R ■ (b X (u X b)) dx dt. 

This allows the geometric depletion to be exploited since, as the alignment between u and b in- 
creases, the value \u x b\ is diminished in a scale-dependent manner. Letting 9(u, b) denotes the 
angle between u and b our interpretation of this phenomenon is the implication, 

H> (£) ^ sin (^ fe ))< (I;) 14 - (10) 

Theorem 11. Let {xj}" =1 C B(xq, Rq) be the centers of a {K\,K2)-cover at scale R. Assume the dynam- 
ical alignment assumption, ( flO) , holds. For appropriate values fa and fa, if 

r := max j (jjj^ ^ ( j|) ^ j < mm{(3 u , fa} =: fa 
then, for scales R where rj (3 < R < Rq, we have, 



2K 1 \Re u Rm V J ~ ' ~ V fie #™ 

Proof. The only procedural difference between the present situation and that seen in (J8)) when 
bounding the quantity, 

(V0) ■ b x (u x b) dx dt. 



To accommodate the dynamical alignment assumption we split the spatial integral over regions 
depending on the relationship \u\ < R/Rq. If this is satisfied then, 

J J (V0) • b x (u x b) dx dt < c (Jpj J J |fc| V p_1 dx dt. 
Otherwise we proceed as in (|9]l but with a minor modification to the normalized radius: 

^ f I p (b x(ux b)) dx dt < (^) 3/4 f [ (K| 1/2 <^- 3/4 )(|6| 1/2 )(|^| 1/2 |&| 3/ V 3/4 ) dx dt 



o 



<co(Mflm) 2 (sup||»||5sup||i.|© ] ^^y jT' |K 4 '- 3 ||| « 



From here one proceeds as in the proof of Theorem[8]to establish fa and fa. The details are omitted 
to avoid redundancy. □ 
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6. Individual Cascades. 

We presently identify conditions under which direct cascades of several individual energy ex- 
change mechanisms - including the exchange of energy between the velocity and magnetic fields 
- are evident in the context of the dimensionless equations. Throughout this section all solutions 
are assumed to be suitable weak solutions to the dimensionless 3D MHD system, regular enough 
for the localized energy equality to hold. 

The localized (by function to a ball of radius R/Ro) flux term responsible for the kinetic to kinetic 
energy exchange driven by the advection of the velocity field is, 

F^ R : =- J(u-V)u-(0u)dxdt = ^ J \u\ 2 (u ■ V<£) dx dt. (11) 

Similarly the (localized) magnetic to magnetic energy transfer driven by the advection of the ve- 
locity field is, 

F h XuR :=- J J{u- Vb) ■ (0 b) dxdt= l -J^ J \b\ 2 (u ■ Vcj>) dx dt. (12) 

The fluid pressure flux-type term is, 

Fl. R := - J Jv P -(4>u)dxdt = J^ J p (u ■ V</») dx dt. (13) 

As already mentioned, the transfer of magnetic to kinetic energy is driven by the Lorentz force 
while the stretching of the magnetic field lines is responsible for the transfer of kinetic energy to 
magnetic energy; combined (locally), they yield the following term - the advection of cross-helicity 
by the magnetic field, 



Until now we have only investigated cascades associated with collections of flux-type terms in- 
cluding the term for the flux of the fluid pressure. Due to the unique structure of this term, ad- 
ditional effort is required to establish cascades of combined fluxes excluding the pressure flux. In 
particular, we will need to bound the quantity, 

l 

p (u ■ V(/>) dx dt. 



(i 



For this, we'll need the following estimate, 

U 1 /2|| 3= ( / (| n | 3 /V 3/4 )(|n| 3 / 2 3 / 4 ) ' 



\3/4/ , \1/4n1/3 

< [ [ / \u\ 2 (j)dx / \uf(j) 3 dx 



< C r ||^0 1 /2||V2|| V(tz0 l/ 2) ||l/2_ 
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Proceeding to the quantity of interest and employing the above bound, we have, 




where we have used Holder's inequality, the Gagliardo-Nirenberg inequality, and Young's in- 
equality. Note that the quantity appearing above, 



/ r l \ 8/9 

C P =(sup||n^ 2 || 2 ) 4/9 (^ \WC 1/a \\l%*) . ( 16 ) 

is dimensionless and a priori bounded in virtue of regularity properties of suitable weak solutions 
(c.f. 11201 ). In addition, it contains no gradients and is independent of the particular cover element 
at scale R. 

Regarding the first term on the right hand side of dl5l ), in order to pass from ensemble averages of 
localized quantities to an integral scale quantity we will need a simple consequence of the finite 
form of Jensen's inequality. Specifically, for {aj}" =1 a set of non-negative values, 



1/9 , n x 1/9 



< 



--1 v t=i 



n \ ' n 



Taking an ensemble average of normalized quantities yields, 

<C P (§) 4 (^)^ (17) 
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Bounds for the remaining flux densities are as follows, 

K,r < f 1 IMblH^llPllvM^)!!^ dx dt 



R 



and, 



t«. < C(Rm"/ 3 sup|M|^ + l)^(f )* [ dt (18) 

Henceforth, we will be concerned with the scenario in which the kinetic energy over the integral 
domain and within the prescribed time-scale remains bounded away from zero. A quick comment 
on the (spatial) scaling of the non-dimensional flow is necessary to give a proper interpretation 
of the following results. Because in the dimensionless context u scales like 1 it is reasonable to 
assume e$ > 1. In the dimensional context we would equivalently assume that el > RlT~ 2 . 
This assumption will be necessary in order to establish cascades for fluxes which exclude the fluid 
pressure flux. 

We are now ready to establish the cascade for the direct kinetic to kinetic energy transfer driven 
by the advection of the velocity field. 

Theorem 12. Assume u and b are solutions of 3D MHD possessing sufficient regularity so that the gen- 
eralized energy equality holds. Let {xj}™ =1 C B(xq,Rq) be the centers of a {K\^K-i)-cover at scale R. 
Suppose that the flow is such that e^ > 1. For certain values /3 U and (3b, if 



max \\w '(it) } <min ^'^> =:/3 ' 



then, for scales R where t / (3 < R < Rq, we have, 



2Ki V Re u Rm V ~ V Re u Rm 



Proof. The main estimates have already been established. Based on the generalized energy equal- 
ity we have 

1 5 

1 f 1 /", ,2 1 f 1 /", ,2 

/ I \u\ <pt dx dt — / / \u\ A(p dx dt 

2 Jo J 2 -^ e Jo J 
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By the assumption on the flow we can replace the exponent of 1/9 in the bound ((17)) with 1. Taking 
ensemble averages, we obtain the lower bound, 

<r>* >- ^»-^(| p^i6^^- c ^Hf)^ 

where 

C u = C(C p Re i/z + Re + l + {MRmf sup | \u\ I2 sup 1 16| 1 1) , 

t t 

and, 

C h = C{Rm A '' i sup||tt||| + Rm + l). 

This is sufficient to establish values for (3 U and /3& (containing no gradients) and conclude in the 
standard fashion. □ 

At this point, noting that F b and F ub both satisfy (|9]), we have already demonstrated all the steps 
involved in proving the existence of cascades for the densities F b and F ub . We consequently omit 
the proofs. The existence theorem for the cascades of the associated energy transfers follows. 

Theorem 13. Let {xj}™ =1 C B(xq,Rq) be the centers of a (Ki,K 2 )-cover at scale R. Suppose that the 
flow is such that eft > 1. For certain values f3 u and /3b, if 

3 « \ V 4 /J> \ V 4 " 



r := max 



{(I) '(I) }<— 



then for scales R where t/(3 < R < Ro,we have, 



K E U o + 4-e£\ < (F b , F ub )n < 2K 2 (±-E% + 



2K 1 \Re u Rm "J- x ' ~ 'V^ e ^ m 

7. Locality of the flux 

According to turbulence phenomenology in the purely hydrodynamical setting, the average en- 
ergy flux at scale R is supposed to be well-correlated only with the average fluxes at nearby scales 
(throughout the inertial range). This phenomenon has been confirmed in lTT6l [8l [TUl [TT^ . 

In the plasma setting, the question of locality has been somewhat controversial. Recently, Aluie 
and Eyink (l| produced an argument in favor of locality of the total energy and cross-helicity 
fluxes. Numerical work also supports locality (c.f. 115] for the case of decaying turbulence). 

In the present setting, locality of the flux is a simple consequence of existence of the corresponding 
nearly-constant turbulent cascade per unit mass. 

We illustrate this on the example of the kinetic energy flux (transported by the velocity). 
Denote the time-averaged local fluxes associated to the cover element B(xi, R) by ^ Xi ,R, 



$x i ,R = ^J o J ^\u\ 2 (u-V&)dx, 



and the time-averaged local fluxes associated to the cover element B(xi,R), per unit mass, by 

(■T 



1/1 / 1 

*x t ,R=ji J o & J ^H 2 (u-V^i)dx. 
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Then, the (time and ensemble) averaged flux is given by 



The following manifestation of locality follows directly from Theorem IT2l 

Let R and r be two scales within the inertial range delineated in the theorem. Then 




In particular, if r = 2 k R for some integer k, 



i.e., along the dyadic scale, the locality propagates exponentially. 




8. A Scenario Exhibiting Predominant m-to-6 Energy Transfer. 



In |IT3l a dynamic estimate is given on the vortex-stretching term (in the vorticity formulation of 
3D NSE) - across a range of scales - using the ensemble averaging process we have illustrated 
above. The purpose was to present a mathematical evidence of the creation and persistence of 
integral scale length vortex filaments by establishing positivity of the ensemble-averaged vortex 
stretching term across a range of scales extending to the integral scale. 

In the induction equation for the magnetic field, the nonlinear term (b ■ V)tt is responsible for the 
stretching of magnetic field lines. Positivity of (b ■ V)u ■ {tfi b) indicates the magnetic field line is 
being elongated, a phenomenon which corresponds to a transfer of energy from the velocity field 
to the magnetic field (negativity would reflect a diminution of the field line and a local transport of 
energy from the magnetic field to the fluid flow). Consequently, to conclude that the predominant 
energy exchange between the velocity and magnetic fields is from the velocity field to the magnetic 
field across a range of physical scales, it will be sufficient to establish (in an appropriate statistical 
sense) the positivity of (b ■ V)u ■ (cf> b) across these scales. Before proceeding to this task we remark 
that recent work (c.f. Ifl5| ) indicates that imbalanced exchanges are common in certain forced and 
decaying turbulent regimes. 

We label the space-time localized quantity of interest as, 



Theorem 14. Let {xi}f =1 c B(xq, Rq) be the centers of a {K\,K2)-cover at scale R. For a certain value, 




if 




then for scales R where r / f3 < R < Rq, we have, 



1 



E b < {V Xi , R ) R < 2K 2 E b . 



2K X 



20 Z. BRADSHAW AND Z. GRUJIC 

Proof. Starting with the dimensionless induction equation (with Rm = S = 1) and assuming b, is 
regular it is routine to obtain, 

(b ■ V)u ■ {b<t> XiiR ) dxdt = J J \Vb\ 2 <p Xi , R dxdt- J ^\b\ 2 (d t </) XiiR - A^r) dx dt 



+ lJ o J \b\ 2 {u-y^ R )dxdt. 



Note that for a refined cut-off function <fi the bounding process evident in the derivation of lO can 
be modified to yeild, 



Consequently, after taking ensemble averages and recalling the last two estimates in ©, 

(V x „r)r < K 2 El + K 2 ( sup | \u\ |t + 2) (^j) " e b 



and, 



2 



(V Xi , R ) R > K 2 E° - _( sup | |u||! + 2) I -^j eg. 
Selecting an appropriate value for j3 allows us to conclude in the standard fashion. □ 



9. A Note on 3D Incompressible Non-Decaying Fluid Turbulence. 

Existence and locality of the kinetic energy cascade in 3D incompressible decaying (no external 
force, global kinetic energy decaying to zero) fluid turbulence in physical scales has been recently 
obtained in llTOllTTIl . 

This was based on existence and locality of the cascade featuring ensemble averaged fluxes of 
time-averaged, cover elements-localized total energy (kinetic plus potential), 

l f T f (W 2 



J / (^2~ + P ) ( U ' V ^ dx dL 



The total energy cascade is universal, i.e., independent of whether global (or integral) kinetic 
energy is decaying or not; the assumption of decaying turbulence was used to obtain a bona fide 
kinetic energy cascade, as postulated in the turbulence phenomenology. 

The estimates on the pressure from Section 6 allow us to now consider the opposite scenario, i.e., 
the scenario in which the kinetic energy associated with the integral domain is bounded away 
from zero - the case of non-decaying turbulence. 

Henceforth, we consider suitable weak solutions to the 3D, dimensionless Navier-Stokes equa- 
tions where, for convenience, Re = 1. We further assume the solutions possess sufficient regular- 
ity to ensure no loss of energy due to potential singularities is evident - i.e., the generalized energy 
equality holds. From this equality we obtain a formula for the flux (localized to the support of the 
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refined cut-off function 0) of interest: 

pi r I |2 pi p pi p I |2 

F< t> = J J ^-(« • V0) dxdt = J J \Vu\ 2 (j) dxdt- J j ^- (<9 t + A0) cte (it 

p(u • V(/>) cfcc (it. 

To proceed note that the bounds from the MHD context on the energy-level localized quantities on 
the right hand side of the above equation all hold. For convenience, recall that these were, 

i 



llo I \ u ^ tdxdt - c °(j[) J j H 2 <P 4p ~ 3 dxdt, 

I r 1 r / p \ 4 r 1 r 

II I' |2a^^^+ ^ „J ^0 1 / / L.|2 jAp-Z 



f \u\ z A(l)dx dt < c l-^ ) j j \u\ z ^ p ~ 6 dxdt, 



f J P (u ■ V*) (ix A < jf* ||^- 3 ||i (ii) 1/9 + C0^V jf 1 ||n^- 3 ||i A 

+ i^ 1 ||(V^ 1 /2||2^ 

where C p is as in (fl6l) . The ensemble average of the sublinear term in the above pressure bound 
is still treated as in (FT)) . That is, upon ensemble averaging localized (by refined cut-off functions 
associated to elements of a (K±, K 2 ) -cover of B(xq,Rq)) dimensionless space-time averages, we 
have, 



n 



where in the last step we imposed an assumption on the flow motivated by the scaling properties 
of solutions to the dimensionless 3D NSE, namely that eo > 1. Ensemble averaging then gives the 
following interpolative bounds, 

1 Eq - C U K 2 (^] e < (F) R < -K 2 E Q + C U K 2 ( e , 



4Ki V-R/ — '"—4 \R 

where, 

C u = C p + C + 2co- 
The correction to the inertial range is now specified to be, 

4 



Provided then that, 



2C U K X K 2 



H,qJ Kq 



our interpolation becomes, for an appropriate constant K, 

^E < {F) R < 2KE . 

Remark 15. The inertial range-correction parameter (3 depends (via C p ) on global a priori bounded quan- 
tities featuring no dependence on Vu (Lf^L 2 and Lf L 3 . norms of u). Consequently, the inertial range 
condition is still a requirement that the gradient of the velocity is large (space-time averaged, over the in- 
tegral domain) with respect to the velocity itself; this will hold provided the integral domain is a region 
exhibiting high spatial complexity of the flow. 
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We summarize our finding in the following theorem. 

Theorem 16. Given the definitions and labelings above, that the flow satisfies eo > 1, and provided R lies 
within the inertial range, i.e., 



we have, 

7^E < (F) R < 2KE . 

Remark 17. The above exposed approach opens up an avenue for the study of forced turbulence in phys- 
ical scales of 3D incompressible flows. This is the topic of the current research project |5|. 
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